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ABSTRACT 

Young stars typically have strong magnetic fields, so that the magnetospheres 
of newly formed close binaries can interact, dissipate energy, and produce syn- 
chrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 
51 day orbit, displays such a signature, with peak emission taking place near peri- 
astron. This paper proposes that the observed emission arises from the change in 
energy stored in the composite magnetic field of the system. We model the fields 
using the leading order (dipole) components and show that this picture is consis- 
tent with current observations. In this model, the observed radiation accounts for 
a fraction of the available energy of interaction between the magnetic fields from 
the two stars. Assuming antisymmetry, we compute the interaction energy E\-at 
as a function of the stellar radii, the stellar magnetic field strengths, the binary 
semi-major axis, and orbital eccentricity, all of which can be measured indepen- 
dently of the synchrotron radiation. The variability in time and energetics of 
the synchrotron radiation depend on the details of the annihilation of magnetic 
fields through reconnection events, which generate electric fields that accelerate 
charged particles, and how those charged particles, especially fast electrons, are 
removed from the interaction region. However, the major qualitative features are 
well described by the background changes in the global magnetic configuration 
driven by the orbital motion. The theory can be tested by observing a collection 
of pre-main-sequence binary systems. 

Subject headings: magnetohydrodynamics (MHD) — binaries: general — stars: 
pre-main sequence — stars: magnetic field 
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Introduction 



Young stars often ha ve strong magneti c fields on their surfaces, with typical field 
strengths 5* ~ 1 — 2 kG fjJohns-KruU II2009I ). When binary star systems have sufficiently 
close separation, the magnetic fields can interact over the course of the orbit. For binaries 
with circular orbits, variability of the radio emission should not correlate with orbital phase. 
For eccentric binaries, however, magnetic interactions reach a maximum near periastron and 
may provide a means for dissipating energy, a portion of which can be emitted as synchrotron 
radiation. 

The T Tauri binary V773 Tau A is an observational example of such a n interacting 
syste m. The stars have masses Mi = 1.54±O.14M0 and M2 = 1.33±O.O97M0 feoden et al 
20071 ). The orbit is observed to have period Porb = 51.1 ± 0.02 days and eccentricity e = 
0.27 ± 0.01. The estimated stellar radii are Ri = 2.22 ± O.2Oi?0 for the primary and R2 = 
1.74±O.19i?0 for the secondary. This binary system is actually the "inner" binary of a more 
extended quadruple system; for purposes of this paper, however, the other two components 
are too distant to affect the dynamics and will be ignored. With the above orbital parameters, 
the semi-major axis = 0.38 AU. The separation 2L between the two stars (twice the 
interaction distance L) equals ao(l + e) ^ 0.48 AU at apoastron and Onfl — e ) ~ 0.28 AU 
at periastron. The radio flaring activity (IMassi et al.ll2006t iTorres et al. 1 1201 if ) varies from 
a few mJy near apoastron, where 2L ^ 50/2*, to a few hundred mJy near periastron, where 
2L ^ 30-R*, if we adopt a mean stellar radius P* = 2i20. 

The scenario that we envision for this source can be described as follows: Each member 
of the binary pair has a strong magnetic field on its surface. Young stars are often observed 
to have magnetic field contributio ns from several mult ipoles, with the octupole component 
being important near the surface (IGregory et al.ll2010l ): however, the dipole field dominates 
at the distance of the interaction region. Over the binary orbit, the distance between the 
stars varies in cyclical fashion. For aligned dipole fields, the magnetic fields become squeezed 
together as the stars become closer, and the magnetic field plays the role of a spring-like 
restoring force; the gravitational force is much larger and little energy is dissipated. For 
anti-aligned dipoles, however, the field lines can connect one star to the other, as depicted in 
Figuredl The field lines that originate at high latitudes connect one star to the other, whereas 
the field lines that originate at low latitudes connect the star back to itself through loops that 
cross the equatorial plane on the sides in the opposite direction from the other star. Note 
that these field lines are much like that of a dipole configuration, but are distorted by the 
magnetic presence of the companion. When the stars are closer, more field lines that start on 
one stellar surface end on the other star; the total energy stored in the composite magnetic 
field configuration is thus smaller at periastron than for larger separations. This trend is 
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illustrated by Figure |2l which shows 10 particular field lines when the system is at periastron 
(top panel) and apoastron (bottom panel). Since the stored magnetic energy changes over 
the orbit, it can provide power for the observed synchrotron emission. Note that energy must 
be given back to the magnetic field configuration as the stars move apart — work must be 
done on the system. Besides dynamo action in the stellar interiors, orbital motion provides 
the ultimate energy source for pumping up the fields after each dissipation episode. Even if 
the energy dissipated in heat/radiation comes entirely from the orbit, however, the implied 
changes in the orbital elements are too small to observe. Although this paper specializes 
to the case of anti-aligned dipoles, we note that reconnection is possible for more general 
magnetic field configurations (so that the analysis presented below is qualitatively valid). 

Previous models of the radio emission from the V773 Tau A binary have been proposed. 
In particular, the emission could ar ise from the interaction of magnetic helmet streamers 



emanating from the stellar surfaces (IMassi et al.ll2006l |2008| ). Our model is similar in spirit 



to this previous work, except in our picture, the helmet streamers are created as a result 
of the interaction rather than pre-exist as independent structures that happen to "bump" 
when the stars pass each other at periastron. The scenario explored herein assumes that the 
energy comes from the large-scale field structure (dipoles), whose interaction characteristics 
are then computed rather than imposed. 



2. Magnetic Interaction Model 
2.1. Formulation 

The two stars of V773 Tau A lack strong infrared excesses, consistent with their clas- 
sification as weak- line T Tauri stars. As a first treatment, we thus assume that the space 
between the stars contains no d y namically signifi cant mass from circumstellar disks or their 
associated winds ( iKonigll Il991t IShu et al.l Il994j ). Instead, we suppose that any tenuous 
plasma that exists above the stellar atmospheres has negligible inertia relative to the elec- 
tromagnetic forces. With the magnetic fields dominant over both matter and gravity, they 
must satisfy the force free- condition that any currents present in the system fiow parallel to 
the lines of force, 

VxB = aB, (1) 

where a is a function of spatial position. Nonzero values of a correspond to non-vanishing 
electric current. The divergence of the above equation, and the condition of no magnetic 
monopoles V ■ B = 0, imply that a satisfies the subsidiary condition B ■ Va = 0, i.e., a is 
constant along field lines. 
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Fig. 1. — Magnetic field lines for an interacting binary system with separation intermediate 
between periastron and apoastron (in the y — plane). Lengths are given in units of i?*. 
The field fines are spaced at uniform intervals on the stellar surfaces. The stars have the 
same size and surface field strengths B^, and have anti- aligned dipole fields. The dashed 
curves show the field lines that mark the boundary between those that close on the star where 
they originate and those that extend to the other star. These limiting field lines effectively 
loop to the other star via z = ±oo. (Only field lines within the plane containing the line of 
centers of the two stars are shown.) 
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Fig. 2. — Magnetic field lines for an interacting binary system at periastron (top) and 
apoastron (bottom) in the y = plane (the plane containing the line of centers of the two 
stars) . Ten field hnes are shown (sohd curves) , where the fields lines start at the same polar 
angles on the stellar surfaces for both cases. The stars have the same size i?* and surface 
field strengths S*, and have anti-ahgned dipole fields. Lengths are given in units of R*. The 
figure shows that a larger fraction of the field lines connect the two stars when they are closer 
together (compare top and bottom panels). The dashed field line (the separatrix) starts at 
different angles on the two stars when the system is at apoastron and periastron. This field 
line has an X-point, where the field strength B vanishes. As the stars press together, loops 
below the X-point reconnect to become the fields above the X-point. 
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Suppose current does not flow to/from infinity, i.e., the system neither gains nor loses 
net electric charge over time. Then a = on open field lines and hence the open field lines 
must satisfy the vacuum-field equations. The closed field lines have two topological types: [1] 
The field lines emerge from the surface of a star and submerge beneath the surface elsewhere 
on the same star. [2] The field lines emanate from the surface of one star and terminate 
beneath the surface of the second star. In either case, in the presence of finite resistivity, 
any currents will dissipate with time and will asymptotically approach zero unless dynamo 
action restores the differential stresses on electrons and ions that led to the currents in the 
first place. Because current dissipation involves the transformation of field energy into heat 
and/or radiation (e.g., synchrotron emission if non-thermal processes of accelerating charges 
are involved), vacuum field configurations represent a lower state of magnetic energy than 
their force-free counterparts. The important point is that current dissipation and magnetic- 
field reconnection tend to enforce a — )> 0, i.e., the closed field lines approach a vacuum-field 
configuration. For the sake of completeness we note that a pure vacuum field configuration 
may not be fully accessible, since finite resistivity only pushes the system towards the lowest 
energy state consistent with its helicity. Here we assume that the helicity is small and 
consider vacuum-field configurations for both open and closed field lines defined by the 
equations 

V X B = and V • B = 0. (2) 

We are thus computing the minimum energy states that would be present in the magnetically 
interacting binary system if current dissipation occurs on a timescale rapid compared to the 
orbital period. This philosophy informs the rest of our analysis until we get to §2.4. 



2.2. Energy of External Magnetic Fields 

As a starting approximation, we take all field lines to be closed, so that the magnetic 
field strength vanishes at infinity. In this case, the most general field configuration associated 
with either star is an exterior multipole expansion. A magnetic multipole of order ^ > 1 
decreases with distance = |r — x^ | from the stellar center as '^~^^^^'^\ where x^- is the stellar 
position (j = 1, 2). In the V773 Tau A system, magnetic interactions occur at distances that 
are large compared to the stellar radii. To leading order, we can ignore all multipoles higher 
than the dipole component £ = 1. The vacuum field configuration of the two stars can then 
be written 

B = Bi + B2, (3) 
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where both stellar fields are pure dipoles with moments nij — rrijZ, 



'■jj'-j 

j 'J '3 



^ = ^ (^2 COS ^ fj + sin e , (4) 



where rj = r — Xj and rj = \rj\. Note that this form for the dipole field and the corresponding 
definitions of the magnetic moments follow Jackson (1962); in particular, rrij = BjR^, where 
Bj is the equatorial magnetic field strength on the stellar surface and Rj is the stellar radius 
(fori = 1,2). 

Let Eb be the energy in the magnetic field exterior to the two stars. If V denotes the 
volume of all space excluding the two spheres instantaneously occupied by the stars, Eb is 
given by 

Eb= f —d'r = f ^d'r + f + f ^^d'r. (5) 

7y87r JySn Jy Sn Jy An ^> 

For each star, let Vj denote the entire space excluding the volume occupied by star j at its 
instantaneous position. The total magnetic energy Eb can be separated, 

Eb — -E'self + -E'int, (6) 

where £^seif represents the magnetic field energy that applies for isolated stars, 

£^self = / + / ^d\ (7) 

and where £^int represents the difference in magnetic field energy resulting from interaction 
of the stellar fields, 

E- [ f "^1^3 f (j^ (8) 

Jy AtT Jy* StT Jy* StT 

where the V* are the spherical volumes of the stars. 

Independent of the stellar motion, Esdf is constant provided that internal magnetohy- 
drodynamics can maintain constant magnetic dipole moments nij within each star. This self 
energy Egeii provides a benchmark energy 

-E'self — Ei + E2, (9) 

where 

D2r)3 pI pod R2d3 

Ej^^j ^di. edi [I + 3/.^] = -f^, (10) 

where the dimensionless variables are fi — cos 9 and ^ = f/Rj- For the system parameters 
adopted below, E'seif = 4.1 x 10^^ erg. 
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The interaction energy Ei^t varies over the orbit and contributes to the observed variable 
radio emission, provided that the orbit is eccentric so that the separation 2L changes with 
time. The final two terms in equation ([H]) for E^^t involve the negative field energies of the 
individual stars in the excluded volumes of the other star; we called the sum of these two 
terms Eexci- The first term involves the mixed dot product, which can be positive or negative, 
depending on whether the magnetic fields from the two stars in the interaction region are 
primarily aligned (Bi ■ B2 > 0) or anti-aligned (Bi ■ B2 < 0); we call this term -Emix- The 
maximally negative interaction energy E^^-i = -Eexci + -E'mix arises when the stellar dipoles 
are anti-aligned: In this case, oppositely directed field lines can squeeze together, reconnect, 
and release energy (ultimately to accelerate particles that create radiation). For the aligned 
case, the fields would dissipate at a (much) lower rate because the current density is lower. 

Here we specialize to the antisymmetric case, where the stars have identical radii -Ri = 
i?2 = R* and magnetic dipole moments of identical strength that are parallel and antiparallel 
to z, the unit normal of the orbital plane, i.e., mi = B^^R^z = — m2. For the sake of 
definiteness, we adopt a standard value = 1.5 kG for both stars. Although Johns-KruU 
(2009) quotes an average magnetic field strength in T Tauri stars of 2.5 kG, this value 
includes all multipoles, so it still might be overly-optimistic to attribute B^, = 1.5 kG to 
the dipole component. However, both stars in the V773 Tau A system are K stars with 
deep outer convection zones, and have short rotation periods under 3 days (Welty 1995), 
properties that correlate with stronger surface field strengths. As a result, the V773 stars, 
which may well be pre-main-sequence progenitors of RS CnV stars, could be expected to 
support unusually well-ordered and coherent (although variable) fields. 

We define 2L to be the instantaneous center-to-center distance between the stars. The 
midpoint between the stars defines the origin of the coordinates. For stars of unequal masses, 
the origin is not at the center of mass and the frame of reference is accelerating rather than 
inertial. Moreover, the coordinate axes rotate at an instantaneous angular velocity Qcm as 
the stars orbit about their center of mass. Given the rotation periods under 3 days, the 
spin angular velocities Uj of both stars exceed Qcm even at periastron. Thus, the system 
will probably support a unipola r inductor with a form reminiscent of the lo- Jupiter system 



(iGoldreich &: Lynden-Belll Il969l ). The induced electric field could accelerate fast particles 
that augment the radio emission, and this process would have a periodicity correlated with 
the orbital phase. However, the changing magnetic field configuration, considered here, 
provides a larger induced electric field (see Section 2.4). As a result, we leave the unipolar 
inductor for future work and focus on the simpler quasi-magnetostatic concepts that we 
believe are responsible for the radio synchrotron emission. 

First we estimate the mixed dot-product term -Emix as follows: For identical stars with 
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anti-aligned magnetic dipoles, one can show that the composite magnetic field fines are 
symmetric witfi respect to tfie mid-plane perpendicular to tfie line connecting tfie stellar 
centers; tfiis plane lies a distance L from eitfier star. We assume tfiat tfie interaction energy 
is given by tfie magnetic energy density integrated from tfiis plane outward to infinity, i.e., 
tfie magnetic energy in tfiis volume is lost because tfie fields from tfie two stars are connected. 
Introducing a form factor /(w, y?, z) to take care of tfie geometric vagaries of tfie dot product 
of tfie mixed fields, -Emix can be written 

Emix = 2 / dip dz Wdw— —5- — = /mean 7T7TT ' (H) 



JL Jo 



wfiere we fiave used a cylindrical coordinate system witfi tfie z axis pointing along tfie line 
connecting tfie stars; tfie mean- value tfieorem allows us to take /(izr, ip, z) out of tfie integral 
and replace it witfi an appropriate mean value /mean- 

Next we estimate tfie energy from tfie excluded volume integrals by replacing tfie field 
strengtfi at eacfi point in tfie volume by its value at tfie center of tfie spfiere and tfien 
multiplying tfiis value by tfie volume of tfie spfiere. If we introduce a factor /exci to make 
tfie result exact, two sucfi integrals yield 

192^/ ■ 

We expect /exci ~ 1 wfien 2L ^ R^,, wfiicfi fiolds for tfie V773 Tau A system. However, tfie 
excluded volume energy -Eexci is smaller in magnitude tfian tfie contribution -Emix by a factor 
of ~ 8(L/i?*)^ ~ 10^ (for V773 Tau A system parameters). Tfie excluded volume energy 
can tfius be neglected to leading order. 



Using tfie above estimates, tfie interaction energy is given numerically by 

2/ n \6/ T \-3 

(13) 



Eint ~ /mean (S.S X lO^^Crg) 



* \^ f ^* 



5* 

ISOOCy yiAxlO^^cmJ V2 x lO^^cm 



Since /mean is positivc for tfie parallel case and negative for tfie anti-parallel case, we expect 
significant release of magnetic energy (yielding enfianced syncfirotron radiation) only wfien 
tfie magnetic dipole moments of tfie two stars are (mostly) anti-aligned. 

Next we determine tfie magnetic interaction energy by numerically evaluating tfie inte- 
gral for tfie mixed dot-product term in equation (E]). Figure [3] plots tfie resulting dimension- 
less interaction energy —Ei^tL^/{B'^Rl) ^ — fraea.nl '^'^ versus tfie fialf-distance L between tfie 
stars. Since tfiis quantity fias little variation witfi tfie binary separation 2L, tfie interaction 
energy scales as E'int oc L^^ to good approximation (consistent witfi tfie estimate from equa- 
tion [H]). Furtfier, tfie inferred value of /mean is close to unity. As a result, tfie maximum 
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energy released over the orbital period of the V773 Tau A binary is 



lapo 
'int 




6.1 X 10^^ erg. 



(14) 



Notice that the energy difference of i^mt between apoastron and periastron — that available to 
convert into radiation — is only a small fraction of E^eM (about one part in 10^). Nevertheless, 
this energy difference is sufficient to explain the observed synchrotron emission. 

As energy stored in the magnetic field configuration varies over the binary orbit, mag- 
netic field lines must change their form. They become increasingly "connected" as the stars 
become closer (see Figure [2]) . The field lines that begin on one stellar surface and end on the 
other stellar surface generally originate near the stellar poles. The stars thus have a "polar 
cap" region where shared field lines begin. As the stars move farther apart (closer together), 
the size of the polar cap decreases (increases). Figure H] illustrates this behavior by show- 
ing the projected polar cap for varying stellar separations. The boundaries of these polar 
regions are obtained by numerically integrating the magnetic field lines outward from one 
stellar surface to determine if the field line intersects the other stellar surface or returns to 
the original one. Note that the polar caps are not symmetric with respect to the geometrical 
pole of the star (that defined by the z-axis), but rather are "tilted" toward the other star 
(which lies far to the right in the figure). More specifically, for each value of y on the polar 
cap contour, there are two values of x. The field line that attaches to the larger algebraic 
value of X connects to the other star through the analogue of the X-point (drawn in Figures 
[Hand [2] for y = 0) in the equatorial plane (also the plane of the binary orbit). The field line 
that attaches to the smaller algebraic value of x connects to the other star through infinity 
(the limiting field line on the back side of the star in Figure [1] for ?/ = 0, i.e., the dashed 
curve that arches higher and higher in z). Finally, we note that Figure H] shows that the solid 
angle A^. of the polar cap is given by the approximate expression A^. ^ 2R^/L (see Appendix 
|A]for a derivation). 



The above discussion estimates the magnetic interaction energy available to convert 
into synchrotron radiation as a function of stellar separation. The corresponding power 
Pmag available for particle acceleration (and radiation) is the negative time-rate of change of 
the magnetic energy and can be written 



2.3. Power from Magnetic Interactions 



dEr 



'int 




(15) 



mag 



dt 



0.1 



0.08 



CO 

* 

, 0.06 - 



.= 0.04 - 



CO 



0.02 - 



10 



20 



30 



Fig. 3. — Dimensionless interaction energy versus stellar half-separation, where L is measured 
in stellar radii and the energy is scaled by Li^. The stars have the same radii and surface 
field strengths, and have anti-ahgned dipole magnetic fields. The scaled interaction energy 
L^Eint varies slowly with L, so that Ei^t ~ L~^. 
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Fig. 4. — Polar cap region for one star of an interacting binary for varying separations. 
The companion star is to the right of the figure. The polar cap is the area of the star for 
which magnetic field lines that originate on one stellar surface terminate on the other stellar 
surface. The contours shown here are projections in the x — y plane. There are actually two 
projections, one for field lines above z = Q (the orbital plane) and one for below z = Q that 
lie on top of each other in projection. Polar cap regions are shown for stellar separations 
ILjR^ — 10 (dashes), 30 (sohd), and 50 (dots). The heavy circle depicts the stellar equator. 
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where r=2L is the instantaneous distance between the stars. The stars execute a Keplerian 
orbit, 

= (16) 

l + £:cos6' ^ ^ 



with semi-major axis Oq, eccentricity e, and orbital angle 6 (IMurray fc Dermottlll999l ). Peri- 
astron occurs when cos^=l. The orbital angular speed follows Kepler's second law, 6 = J/r"^, 
where the specific angular momentum J = ao^^{l — e^)^^^, and Vt = 27r/Porb is the mean 
angular velocity of the orbit (Porb ~ 51 days). Collecting the above results, the power 
becomes 

_ nBlRl f l+ecosO Y esine 

/mag -/mean \^ ^ _ ^2 J (1_ ^2)1/2- ^'') 

The positive sign in equation f lT7|) indicates that magnetic energy is released as the stars 
become closer together. 

Figure |5] shows the magnitude of the total power as a function of time, which is measured 
in days since apoastron. The plot presents the dimensionless power |Pmag/'^o|, where we have 
defined the fiducial power scale 

= M^i^ ~ 1.3 X 10^9 erg s-\ (18) 

Note that the power vanishes at periastron when sin 6 = 0, whereas peak power occurs ~ 4 
days earlier. For simplicity, in the figure we assumed | /mean I = 1- 

The synchrotron emission is not necessarily fully correlated with the magnetic dissipa- 
tion: Although the deformation of the magnetic field is slow and steady, with timescales 
characteristic of periastron passage (i.e., days; see Fig. |5]), the energy release from the mag- 
netic fields and the resulting radio fiares can take place on significantly shorter timescales. 
For example, the Sun twists up its surface fields on a natural timescale given by the differ- 
ential rotation between its equator and pole (many weeks); when outbursts arise, however, 
they are much more rapid, lasting only minutes for impulsive fiares and hours for coronal 
mass ejections. This phenomenon is much like the buckling of a metal plate stressed along 
its edges: the accumulation of stresses can be slow, but the relief of those stresses — after 
they become supercritical — can occur on catastrophically fast timescales. 



2.4. Acceleration of Particles to Relativistic Energies 

The loss of magnetic field energy does not directly result in radiation because the orbital 
frequencies, or even peak rates of reconnection, are too slow to result in meaningful electro- 
magnetic radiation. Instead, when magnetic energies are changing because of annihilation 
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Fig. 5. — Power dissipated by the magnetic field configuration of the binary system over one- 
half orbit. Time is measured in days since apoastron. The dimensionless power I'Pmag/'^ol is 
plotted on the ordinate, where the fiducial power scale Pq is defined by Vq = B'^R^Q/ao^ ~ 
1.3 X 10^^ erg s~^ (see equation [18]). Energy is dissipated over the incoming portion of the 
orbit. 



- 15 - 



of magnetic field B, electric fields E arise that are given by Faraday's law of induction: 

VXE=-If, (19) 

with c equal to the speed of light in a vacuum. Although the resulting electric field is easy 
to compute as a function of space and time in the vacuum-field approximation, we do not 
record the results here. Wc merely comment that this electric field E does have a component 
parallel to the magnetic field B, and that B ■ E along magnetic field lines changes sign over 
the orbit cycle, so that the electric currents that arise in response to the electromagnetic 
fields are alternating in character and do not involve a permanent charging up of either star. 

The resulting electric field E that is not strictly perpendicular to the magnetic field 
B then accelerates particles, in particular electrons, and it is these fast electrons, if they 
achieve relativistic energies, that emit synchrotron radiation in the ambient magnetic field. 
Some of the fast electrons may have mirror points below the surface of one of the stars in 
the system. In that case, electrons can spiral to the footpoints of the magnetic field and 
collide with ions in the dense atmosphere of the polar cap of that star. Such collisions 
produce brehmsstrahlung radiation, which provides another channel through which the loss 
of magnetic energy can appear as electromagnetic radiation. 



2.5. The Source Region for Fast Electrons 

In a fully ionized plasma, such as may exist in the upper chromospheric regions of the 
two weak-lined T Tauri stars, the effective frequency at which electrons are slowed down 
by Coulomb collisions with other charged particles decreases with increasing energy ^ymeC^, 
where nie is the rest mass of the electron and 7 = (1 — v^/c^)'^^^ is its Lorentz factor, 
and where v is the electron velocity. As a consequence, in the presence of a non-vanishing 
component E\\ of the mean electric field E that is parallel to the mean magnetic field B (i.e., 
£"11 = E ■ 6 where b = B/|B|) electrons above a certain threshold energy can run away to 
ever higher energies (Dreicer 1959): 

meC^ (P-) = -eE • V, (20) 



dt 



elec 



where e is the charge of the electron and the important component of v in this context is 
that parallel to B. 

The ratio of the magnitude of the quasi-electrostatic force to the gravitational force 
at the poles of either star of a particle with charge e and mass nie is e\E\\\Rl/GM^me — 
r\E\\/B^\ where F = eB^Rl/GM^rUe ~ 10^°. Thus, when |£;||| is larger than ~ lO'^^B^, 
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which is true of the vacuum field values, the electrostatic acceleration of electrons (which 
drag ions with them) off of the poles of one star (or the other) will more than overcome the 
gravitational attraction of the stars for the plasma. 

Space charge can accumulate and AC current fiow can occur that modify the self- 
consistent fields relative to the vacuum approximation of two (naked) moving magnetic 
dipoles. The synchrotron radiation observed in the system suggests that there are Ue ~ 
2 X 10^ electrons per cm^ in a volume of space Lr" when the system is at periastron L = 15-R*. 
If these electrons all came from two polar caps, each of area SnRl/L supplied at a current 
flow / ~ {IGn Rl/ L)enp*c over a time equal to half an orbital period Porb/2, electron number 
conservation requires /Porb/2 = eUgL^- In other words, the number density of fast electrons 
in the source region of the polar caps must be Ue* rie(L^/87ri?^cPorb) ~ 426 cm~^. The 
magnetic field associated with the current fiow Bac ~ ene*P*/47r ~ 2300 G is comparable to 
the strength of the dipole magnetic field 2S* = 3000 G at the stellar pole. Thus, we expect 
the back reaction of the induced current fiow of fast electrons to interfere by an order unity 
amount with the vacuum approximation. 

An important quantity in this problem is the number density rifast of fast (runaway) 
electrons that develop when the thermal plasma of a local number density Ue in the chro- 
mosphere or corona of one of the stars is subject to acceleration by an electric field E\\ of 
appropriate sign parallel to the polar-cap magnetic field of strength 2P*. Following the treat- 
ment of Tandberg-Hanssen & Emslie (2009), we can write the rate of fast electron production 
in the form 

0.35neZ/J(e) , (21) 



dt 

where e = |/|E£)|, the strength \Eid\ of the Dreicer electric field is given by 

|Ec| = -^lnA, (22) 

and where this treatment is valid in the limit e ^ 1. The parameter Ad = [kT /AT^Uee^Y^"^ 
is the Debye length and the Coulomb logarthim A takes the form InA = 18 -|- ln(3.16 x 
lQ-52-3/2^-i/2^_ Finally, Pc is the electron-ion colhsion-frequency, 

i/c = 3.6 T-^/V In A, (23) 

where all quantities are in cgs units. In the regime e < 1, the function /(e) has a simple form 
(Tandberg-Hanssen & Emshe 2009); however, for the typical parameters of this problem 
(ne ~ 4 X 10^ cm-3 and T ~ 2 x 10'^ K), the Dreicer field |Ed| ~ 2 x 10"^^ statvolt/cm, so 
that e ^ 1. In this regime almost all of the electrons are subject to runaway. 



The resulting fiow of current, which is almost instantaneous on the orbital time scale, will 
charge up both stars (one negatively and one positively), and nothing in between, because 
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the induction electric field satisfies V • E = 0. The electrostatic potential associated with 
the negative charge grows to provide just enough repulsion of electrons (with a similar 
development for the positive charge on the other star in repelling ions) so that the AC 
current automatically shunts back and forth to maintain a quasi-steady state. In other 
words, a large-scale dipole electrostatic field (with one pole on each star) develops so that 
it compensates macroscopically for the induced electric field of the two moving magnetic 
dipoles. Because of the different spatial dependence of the fields of two moving electric 
monopoles and two moving magnetic dipoles, the cancellation of electric fields only holds as 
an average for the polar caps and not throughout space. Acceleration of electric charges in 
between the two stars occurs in response to the superposed magnetic and electric fields. In 
what follows, wc provide order of magnitude estimates by considering only the action of the 
induction electric field. 



2.6. Constraints Set by the Observed Synchrotron Radiation 

To make a numerical estimate for the electron acceleration that occurs in the spaces 
between the two stars after the electrons leave the source and sink regions near the stellar 
poles, we need to specify a "typical" field strength 5int that is undergoing change on a time 
scale ts over a length scale L in the interaction region. As a typical -Bjnti we therefore take 
the average of its value at the X-point (where B = 0) and its value halfway between one of 
the stars and the X-point. On a direct line, this means that the distance from one star is 
L/2 and from the other star, 3L/2. For anti-aligned dipoles, the magnetic field strength at 
the latter point equals 



Bi/2 — B^R\ 



L/2 J \2>L/2 



ff (^ 1 . (24) 



For our standard parameters, ~ 3.96 G. An average of this value and the X-point value 
{Bx — 0) produces the estimate 

Ant = 2 G. (25) 
The magnitude of the electric field is given by 

1^' ^ (^) ^-'^ ^^^^ 

where r is the time scale on which the magnetic field changes. This time scale is bounded 
from above by the orbit time Porb ~ 4.4 x 10^ s and is bounded from below by the burst 
time scale ts — 10^ s (from the observations of Massi et al. 2006, 2008). For all possible 
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values of the time scale r, the length scale cr ^ L = 2 x 10^^ cm and hence E ^ Bi^t in cgs 
units. Using a (conservative) value near the upper end of the range, r = 10^ s, we find that 
cr = 3 X 10^® cm and |E| ~ 10~^i?int. Nevertheless, as we shall see, the action of this electric 
field over the long distances available in the system can achieve impressive accelerations of 
charged particles. 

In this work we assume that the electric field E generated through Faraday induction 
(from equation [2B]) is the dominant contribution. In particular, we verify below that |E| is 
larger than the field of a unipolar inductor, where an electric field is produced by the slippage 
of magnetic fields through the stellar atmosphere(s). This latter field strength, denoted here 
as |E|uni, can be estimated through the expression 

where we estimate the speed v = |v| of slippage as 

v^J^-—, (28) 

r 



rot 



where Prot is the rotation period of the star, and where the dimensionless fraction J-" ~ 1 
because the orbital angular speed of the companion star is small compared to rotation speed 
of the star itself. Note that to obtain this estimate, we evaluate the magnetic field strength 
from one star at the position of the other star. Comparing this result to the electric field 
obtained from Faraday induction via equations fl2^ - [26l) . we find the ratio 





E| 




|E 


'|uni 



^^~'~ 20, (29) 



2171 J^R^ T 

where we have used r = 10^ s and J-" = 1. The electric field from Faraday induction (used 
here) is thus safely larger than that of the unipolar inductor, although the latter could be 
included correction in future work. 

The maximum frequency of the synchrotron radiation produced by an electron with 
Lorentz factor 7 spiraling in a magnetic field of strength B (see, e.g., Shu 1991) is given by 

V = '^'^vl where vl = 77—— > (30) 



i.e., ul is the Larmor frequency. For B = B\^t = 2 G, in order to produce = 90 GHz 
radiation, 7 must be at least 127. Such a Lorentz factor corresponds to moderately relativistic 
electrons, and we may approximate v ^ c when v stands alone, such as in equation fl20|) . 
We now consider the length of time that a fast electron on field lines connecting the two 
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stars can be expected to experience an electric field such that E • v < so that it receives 
continuing boosts of energy via equation f l20|) . If the fast electron shuttles back and forth 
between mirror points, the maximum length of time that it spends accelerating rather than 
decelerating is ~ 2L/c. If it is doomed to strike a polar cap before mirroring, the maximum 
time is again ~ 2L/c. Integrating equation (120|) over the time interval ILjc^ we find that 
the maximum boost yields a Lorentz 7 given by 

In fact, such high values of 7 are inconsistent with limits obtained from considering the 
energy drain due to synchrotron losses. The maximum 7 consistent with synchrotron losses 
on a time scale r = is given by 

7max~^^^~ 26,000, (32) 

where gt = Snr'^/S = 6.65 x lO^^^cm^ is the Thomson cross section with re = e^/me(? 
equal to the classical radius of the electron (see, e.g., Shu 1991). The maximum value 
7max ~ 26, 000 is considerably larger than the value 7 ~ 127 required to produce the 90 
GHz synchrotron radiation. Having the maximum value 7max larger than the required value 
of 7 is not necessarily a problem, but our arguments below suggest that the synchrotron 
spectrum does not extend to frequencies much higher than 90 GHz, so that some reduction 
of 7max may be indicated. 

Toward this end, we speculate that the self-consistent electric fields are substantially 
reduced by the collective response of the intervening semi-relativistic plasma, and that the 
energies of fast electrons have to be renewed on each periastron passage because the particles 
escape from the system or are lost to brehmsstrahlung as the mirror points change with 
reconnecting fields every orbit. The relevant decay time for the bulk of the fast electrons in 
the system is then r ~ Porb = 51.1 d ~ 4.4 x 10^ s. The electrons that can carry over from one 
outburst to another have a maximum Lorentz 7max ~ 59, which implies that the electrons 
responsible for the 90 GHz radiation would need re-generation for each orbit cycle. A third 
possibility is that the component of E along magnetic field lines is considerably smaller 
than its absolute value |E|, so that an integral of equation (!20|) produces results appreciably 
smaller than the naive estimate of equation (13T!) . We are in the process of quantitatively 
evaluating the cumulative result of all three effects by doing test orbit integrations and will 
report on the results in a future communication. 

From where then does the time scale Ib = 10^ s arise? The time scale for magnetic 
reconnection is given empirically in the case of solar flares roughly by the formula: 

tB = — , (33) 

qvA 
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where va = -Bint/(47rp)^/^ is the Alfven speed for a plasma with mass density p and q 
is a dimensionless number (the Petschek parameter) whose exact value is of considerable 
theoretical debate but probably has a value between 0.1 and 1. In order to achieve the burst 
time scale ~ 10^ s with L = 2 x 10^^ cm, we require an Alfven speed va = 2000g~^ 
km s~^. With our estimate that Ant = 2 G, we then deduce a typical plasma density of 
p ~ 8.0 X 10~^'^g^ g cm~^. For a fully ionized gas with solar abundance, such a plasma would 
have a number density of electrons given approximately by rie ~ 4.2 x lO^g^ cm~'^, well 
below the value, for g < 1, that would have caused (an unseen) Faraday de-polarization of 
the emitted synchrotron radiation (Phillips et al. 1996). Note that the communication speed 
f A is large compared to the orbital speeds at which the magnetic configuration is stressed; 
this ordering justifies our treatment of the background magnetic evolution as being governed 
by elliptic partial-differential equations (the vacuum-field equations). 

We still need an explanation for the value of the number density Ue ~ 4.2 x lO^g^ cm~^. 
The flux density at a frequency u = 90 GHz during the radio outburst reaches a temporal 
peak (see Figure 2 of Massi et al. 2006) given by 



followed by an exponential decay with mean life 2.3 hr, which we interpret to be essentially 
the reconnection time scale ts = 10^ s, and not the radio-synchrotron decay time scale 
T ~ -Porb = 51.1 d. We compute the total emitted synchrotron power from the product of 
the peak Si, and t_B = 10^ s integrated over all frequencies. 

To perform the integration over u we have to specify the energy distribution of the 
relativistic electrons. In our situation, the distribution is established, not by a coUisionless 
shock, but by the acceleration provided by magnetic reconnection that varies with the binary 
orbit balanced by losses through particle removal from the polar caps (and the X-region; see 
below). For this problem, we lack an a priori theory for the number density n(7) d'j of fast 
electrons with Lorentz gamma between 7 and 7 + c/7. In Appendix |Bl we discuss the case 
of a general power-law distribution. 



For the sake of definiteness, we start by specializing to the case with p = 2.4 (the value 
appropriate for coUisionless shocks). For p = 2.4, the integral of n{'~f) d'y from 7 = 1 to 
7 = 7max is insensitive to the cutoff value ■jrna.x as long as 7max ^ 1 (e.g., if 7max = 127). The 
total number density of fast electrons in such cases is approximately no/1.4. 

The quantity no yields the volume emissivity of synchrotron radiation as 



S„ = 400 mJy = 4 x 10"^^ erg s"^ cm'^ Hz 



-1 



(34) 



n(7) d'y = 77.07 ^ . 



(35) 




(36) 
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where r^. is again the classical radius of the electron, v± = eB±/2'n-'meC, and B± is the 
component of the magnetic field perpendicular to the line of sight (i.e., in the plane of the 
sky for angularly resolved observations). The dimensionless number Ci ^4.1 for our choice 
p = 2.4, and is defined in general by equation flB4p in Appendix [Bl For optically thin 
synchrotron emission, the fiux density at a distance d from the source is given by 



where we have used equation (18.20) of Shu (1991) [which is missing a factor of Rs = L], 
and where L is the characteristic "radius" of the source. 

With S,y (X z/^°'^ in equation ( l37|) . the implied integrated power from zero frequency to 
z/ is given by 



where d = 133 pc is the distance to V773 Tau (Torres et al. 2011). If we put in the 
observed value Si, = 400 mJy = 4 x 10~^^ erg cm^^ s^^ Hz~^ at z/ = 90 GHz, we find that 
Vu = 2.5 X lO^'' erg s~^ for the synchrotron emission from zero frequency to z/ = 90 GHz. 
Over a burst time of = lO'' s, the released energy at radio frequencies up to 90 GHz equals 
2.5 X 10^^ erg, which is nominally 40% of the total magnetic energy available in our model 
in going from apoastron to periastron, E^^° — Ef^^^ = 6.1 x lO^^erg (refer to eq. [I4]). 

The above considerations suggest that if = 1500 G is a good estimate for the stellar 
surface fields, the conversion of magnetic energy into synchrotron radiation must be efficient, 
and little synchrotron radiation is emitted at frequencies higher than 90 GHz. Because the 
collisional losses to the thermal plasma are expected to be negligible (see the discussion 
in Appendix A of Massi et al. 2006), we believe the acceleration of runaway electrons to 
be highly efficient, not only in the sense that a high proportion of eligible plasma can be 
driven to relativistic energies, but also to the extent that there is a limited escape of such 
accelerated electrons from the system before they are overtaken by synchrotron losses. In 
other words, as long as the cosmic rays generated within closed field lines do not open such 
field lines, almost the entire energy available from the change of the magnetic configuration 
will be used to accelerate the particles responsible for the deduced radio emission. Moreover, 
these particles essentially lose all their energy by synchrotron emission before the next cycle 
of field dissipation, particle acceleration, and synchrotron emission; as a result, almost all 
the energy input from magnetic field dissipation in the V773 Tau A system is converted 
into radio synchrotron emission. There is, however, a supra-thermal population of electrons 
that carries over from one cycle to another that serves as the seed population (the "storage 
ring") for the next round of synchrotron acceleration. In fact, we argue below that the 




(37) 




(38) 
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supra-thermal population may be the predominant fraction of the plasma density that we 
estimated using the reconnection time argument. 

For Bj_ = 2 G and z//z/_l = (127)^, we compute from equation flTTj) and 5*;^ = 400 mJy at 
u = 90 GHz that no/1.4 = 1.44 x 10^ cm^. Within the uncertainties of the calculation, the 
number density nfast of fast electrons needed to account for the observed radio emission is 
then comparable to the total electron density Ue = 4.2 x lO^g^ cm~'^ deduced to be present on 
the basis of the reconnection time. (The latter calculation was performed non-relativistically, 
which is justified for the ions that provide the bulk of the mass density. The number density 
of the electrons, which may be relativistic, is then deduced on the basis of charge neutrality.) 
In other words, the supra-thermal population may be a significant fraction or even all of the 
plasma electrons (if q = 0.19), which provides the sought-for explanation for the value of 
rig. As a bonus, we obtain the attractive suggestion that the supra-thermal electrons may 
themselves be the contributor to the "anomalous resistivity" that triggers the reconnection 
event. 

That such a source of anomalous resistivity is needed can be seen by computing the 
reconnection scale if it occurs by collisional resistivity: iji = (t^^?)^^^ ~ 3 km (T/lO^i^)"^/^, 
where we have used the Spitzer (1978) formula for the resistivity of a fully-ionized cosmic 
plasma: i] = lO^^T"^/^ cm^ s~^. For any plausible value of the temperature T of a thermal 
plasma, the length scale in over which collisional resistivity can induce reconnection is much 
smaller than the scale L over which the spatially-resolved VLBI observations require the 
events to occur (Massi et al. 2006). We conclude that the true resistivity occurs not by 
physical collisions knocking charged particles off field lines, but by their scattering off fluc- 
tuating electromagnetic fields that arise because the laminar collisionless flow (particularly 
the counter-streaming of charges driven by quasi-electrostatic fields from the polar caps of 
the two stars) is unstable to the generation of plasma waves. 

We can now provide a consistency check on our deduced value of no using the following 
mechanical argument: If we follow convention and assume that the momentum distribution 
of cosmic-ray electrons is isotropic, we can define a cosmic-ray electron pressure Per through 
the expression 

Per = - / vpn{-f) d-f ^ -nieC^no / 7"^'^ 6/7 ^ 0.83meC^no, (39) 
3 Jo 3 Jo 

where we have extended the integration to infinity because most of the contribution comes 
from low energies, and where we have made the approximation that the product of speed and 
momentum, vp, for a relativistic electron equals its energy 'jniec'^- Using no/1.4 = 1.44 x 10^ 
cm~^, we find that the electron pressure P^ ~ 0.14 erg cm~^. This value should be compared 
with the magnetic pressure Pj^t/87r ^ 0.16 erg cm~^. The pressure of fast electrons at peak 
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outburst is thus comparable to the typical magnetic field pressure in the interaction region, 
implying that the electrons that escape from beneath the X-point into the region above the 
X-point in Figure |2] can infiate the fields in the latter until they become open. Once every 
orbit cycle, the fast electrons can escape their confinement through such a mechanism; as a 
result, the density no cannot accumulate through repeated accelerations to a value that is 
much in excess of that calculated here. Similar considerations in other contexts may offer an 
explanation for why the assumption of rough equipartition of energies in cosmic-ray electrons 
and magnetic fields holds in many systems of interest in astrophysics. 

For sufficiently low frequencies, synchrotron emission is suppressed due to the Razin- 
Tsytovich effect (e.g., Razin 1960, Tsytovich 1951, Simon 1969, Melrose 1972). The lowest 
possible value of the frequency for which waves can propagate is given by 

= TT^^ - 20^^^Hz ~ 42 g^MHz , (40) 

where up = {uec'^ /TirneY^'^, col = eB/nieC, and where we have used Ue = 4.2 x lO^g^cm"^, 
B ~ 2G, and sin^ ~ 1. This frequency is much lower than the 90 GHz synchrotron emission 
observed in the V773 binary system, so that Razin-Tsytovich suppression does not appear. 

We check now whether the modeled synchrotron emission is optically thin as assumed. 
The mean synchrotron optical depth through the source region is given by 

where the constant C2 ~ 9.9 for p = 2.4, and is defined in general by equation (]B9P in 
Appendix [Bl The right-hand side of equation (HTl) nominally equals 5 x 10~^ at 90 GHz, 
which would indeed make the synchrotron radiation optically thin at the observed radio 
frequency. Radio observations at lower frequencies could be used to measure the spectral 
index and to show the bend in the spectrum as a transition is made to optically thick 
conditions; such measurements would yield valuable diagnostic information on the physical 
conditions in the system. 

Finally, we note that the ratio of inverse Compton losses to synchrotron losses is given by 
u^/ub where u-^ is the energy density of photons and ub = -B^/Stt. Most of the energy density 
in the V773 Tau A system arises from optical photons emitted by the t wo stars; to estimat e 



this contribution, we divide the optical luminosity of the system, 3.93 Lq (IBoden et al.ll2007l ). 
by AtiLF'c to obtain ~ 0.010 erg cm~^. This value should be compared with the magnetic 
counterpart ub ~ 0.16 erg cm~^. Although inverse Compton losses are significantly smaller 
than synchrotron losses, they are not completely negligible. We suggest that it might be 
worthwhile to try to detect optical photons that have been scattered into the hard X-ray 
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regime by the inverse Compton effect involving tfie syncfirotron electrons that are deduced 
for the V773 Tau A system. 

In summary, if our model is to have validity, two conditions must hold: 



• The synchrotron spectrum Si, must not extend to frequencies much higher than the 
observed radio frequency of = 90 GHz. At first sight, this requirement seems to be 
contradicted by the case of the DQ Tau binary, where an X-ray outburst coincided with 
a radio outburst (Getman et al. 2011; see also Salter et al. 2008). However, enhanced 
X-ray emission also occurs in the case of solar radio outbursts, and the interpretation 
there is that the fast particles created by electromagnetic processes in an impulsive 
flare or coronal mass ejection travel down the field lines to the footpoints of a magnetic 
loop. The collisions that the fast electrons have with the solar atmosphere then produce 
optical and X-ray emission by brehmsstrahlung radiation. It would be interesting to 
know whether the X-ray outbursts in DQ Tau are synchrotron events (i.e., with a power- 
law spectrum and a high degree of linear polarization), or they are brehmsstrahlung 
events (i.e., with an optically-thin thermal spectrum and no associated polarization), or 
they have the characteristics of high-energy photons that have been created by inverse 
Compton scattering. 

• Although the change of magnetospheric configuration in the V773 Tau A binary is 
occurring on the orbital time scale of weeks, the actual release of the stored energy 
occurs not on the orbital time scale (i.e., not as depicted in Fig. 4), but on a much 
shorter burst time scale of hours (e.g., Ib ~ lO'' s). Again this property is consistent 
with the solar experience. Although the twisting of the coronal magnetic configuration 
originates in differential motions at the footpoints of the field in a solar atmosphere 
whose rotation rates at the pole and the equator are measured on a time of many weeks, 
the outbursts when they come occur on a time scale of minutes (for impulsive flares) 
or hours (for coronal mass ejections). The build-up of magnetic stresses may occur 
slowly and steadily [without instantaneous dissipation as assumed by the simplifying 
philosophy adopted at the beginning of our discussion (see §1)], but the rehef of such 
stresses can appear suddenly and catastrophically as in solar flare events. 



2.7. Stellar Winds 



If the stars in the system have substantial stellar winds, their mechanical luminosity 
could compete with magnetically produced power from the interaction region. Scaling the 
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mechanical luminosity Lmech of the wind to the outflow rate from the Sun, we find 

If the stars in V773 Tau A had solar-type winds, their mechanical luminosity would be more 
than 200 times smaller than the magnetic power scale of equation ( !T8|) . With a mass- loss 
rate M ^ 10~^^ Mq yr~^ for each star, wind power would be competitive with magnetic 
power, and such mass-loss rates are not beyond the capabilities of weak-line T Tauri stars. 

In settings where stellar winds help shape magnetic field structures in young stars, 
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more slowly than a magnetic dipole, and coUisionless shocks might play a stronger role in 
accelerating pa rticles than magn etic reconnection. However, VLBI mapping for the V773 



Tau A system (IMassi et al.ll2008l ) suggests otherwise; these authors interpret the magnetic 
field configurations as solar-like helmut streamers, where the outflow is not strong enough to 
overwhelm the fields. As a result, we suspect that the mass- loss rate in stellar winds in each 
of these weak-line T Tauri stars is substantially smaller than M ^ lO^^^M© yr~^. Indeed, 
observations suggest that typical mass loss rates from weak-lined T Tauri stars are roughly 
M ~ lO-i^M© yr-i (Guenther & Emerson 1997). 

Although the mechanical luminosity of the wind could compete with the magnetic power 
dissipated in the interaction region, the wind is not strong enough to greatly alter the 
magnetic field structure: The magnetic pressure in the interaction region is given by Pb = 
where we expect Ant ~ 2 G, so that Pb ~ 0.16 dyne cm~^. For comparison, the 
ram pressure of the wind, evaluated at the location of the interaction region, can be written 
in the form Pram = Mf^/(47rr^). For a mass loss rate M = IO'^^Mq yr~^, the ram pressure 
is only Pram ~ 5.6 x 10~® dyne cm~^, nearly 30,000 times smaller than the magnetic pressure. 
Closer to the star, the magnetic pressure dominates by an even larger margin. As a result, 
the magnetic field is strong enough that the wind cannot change the field configuration, and 
hence the use of dipole fields remains valid. 



2.8. Back Reaction 

In this scenario, magnetic energy is released when field lines connect as the stars move 
closer together. After periastron, as the stars separate, energy must be converted back into 
the magnetic fields in order to maintain the same dipole strength. This energy could come 
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from either the stellar orbit or the stellar spins. Using equations f lTTj) and f|T6l) . we can 
estimate the total energy released per orbit 

^ ' , (43) 



3al J (l-e2)3^ 

where Oo is the semi- major axis and we have assumed |/mean| ~ 1- Using parameters from 
equation (fT3|) . we find the change in energy per orbit AE ?s 6 x 10^^ erg (consistent with eq. 
|14]). The available mechanical energy in the orbit is the difference between the energies of 
the eccentric orbit and the circular orbit with the same angular momentum, about 3.5 x 10^^ 
erg. This energy supply can last 6 x 10^° orbits, much longer than the age of the system. 
The back reaction on the stellar orbit is thus negligible. 

Next we consider the stellar spins, which provide another possible source of energy for the 
magnetic interactions. Field lines connecting the two stars will be twisted by the combined 
motion of the orbit and the stellar spins to produce an azimuthal field. After one rotation 
period, these azimuthal field lines must reverse polarity along the line of centers joining the 
stars. This magnetic configuration naturally generates current sheets which in turn trigger 
reconnection events (because the field lines cannot be sheared indefinitely), consistent with 
the picture developed in this paper. Since the stars have periods of approximately Prot ~ 3 
day, both bodies have total spin kinetic energy K^pi^ ~ Jw^, where the moment of inertia 
/ = kM^Rl ~ k{Q X 10^^) g cm^, and where is a dimensionless constant. For example, for 
a polytrope with index n = 1.5 (a good approximation for completely convective pre-main- 
sequence stars), one finds k ^ 0.205. The total energy stored in stellar rotation is then given 
by -fCspin ~ 7 X 10^^ erg, which is comparable to, but somewhat larger than, the available 
energy stored in the orbit. As a result, the spin energy supply can last for billions of orbits 
and the backreaction can (again) be neglected. 



3. Conclusion 

This paper shows that interactions between the magnetic fields of young eccentric bina- 
ries can provide significant power for accelerating fast electrons. The available interaction 
energy is large enough to produce the observed radio synchrotron signature in the well- 
studied system V773 Tau A. The energy stored in the magnetic fields of the system varies 
with orbital phase, with a maximum rate of energy dissipation occurring just before peri- 
astron (§2.3, Fig. E]). This basic model can explain qualitatively the synchrotron radiation 
observed in the T Tauri binary V773 Tau A. However, the details as to how and why the 
radio flares occur with the energetics and time scales that they exhibit involve more elab- 
orate considerations. These considerations include how magnetic field annihilation yields 
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electric fields that accelerate electrons to relativistic energies and how these electrons are 
ultimately lost to the system, through synchrotron radiation, or brehmsstrahlung events, or 
escape from the system. 

In the last regard, we began by assuming that all the field lines in the system close onto 
one star or the other. However, in §2.4 we presented the case that during peak outburst, 
the fast electrons in the V773 Tau A system may have enough pressure to open some of 
the field lines that join one star to the other. Through electrostatic forces, the flow of the 
electrons to the interstellar medium will carry a corresponding number of ions with them. 
The accelerated fast electrons that do not escape from the system will ultimately lose energy, 
either through collisions with ions (mostly at the polar caps of the two stars) or through 
synchrotron radiation. A surprising prediction of this model for the V773 Tau A system is 
that the synchrotron spectrum cuts off at radio frequencies only somewhat higher than the 
observed value of 90 GHz. Another prediction is that the X-ray emission in th e young binary 



DQ T au (with periastron ~ 8-R*), which is correlated with the radio bursts (IGetman et al. 



201ll ). arises not from synchrotron emission, but probably has a brehmsstrahlung origin 
when fast particles associated with the electromagnetic acceleration stream to the footprints 
of the magnetic field and strike the surface of one of the stars in the system. 

In closing, we note that the ideas presented in this paper can be tested. Perhaps the 
best way to study these phenomena further is through observational surveys of synchrotron 
radiation in a large sample of pre- main- sequence binary systems. In this model, the magnetic 
interaction energy E^^t and the idealized dissipated power Pmag are straightforward functions 
of the stellar field strengths, stellar radii, semi-major axis, and orbital eccentricity (eqns. 
[TT] . [13], and [IZ]). The binary orbital elements and stellar properties can be measured 
independently of the synchrotron radiation. While hysteresis will make the synchrotron 
power more sporadic than the idealized dissipation described by equation (ITTl) . the total 
energy release per orbit should be derivable from equation 
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A. Solid Angle of the Magnetic Polar Cap 



As discussed in tlie text, and shown in Figure HI the solid angle Ac of the polar cap is 
given by the approximate empirical formula Ac ~ 2R^,/L. This result can be understood as 
follows: Half of the circumference of the polar cap arises from field lines that connect to a 
straight line of X-points defined hj {z = 0,x = L), with y running from — oo to +00. In 
spherical polar coordinates {r,9,ip), the equation for a magnetic field line of a single dipole 
placed at the center of a star (the origin of the coordinate system) is defined by the equation 



ro sin^ 9, 



(Al) 



where the field line crosses the equatorial plane at r = tq. The line of X-points is given by 
the equation rocosip = L with ip running from —tt/2 to +tt/2. These lines map onto the 
surface of the star r = through the equation 



sin"* = [ ^ ] cos (p 



(A2) 



which describes a semi-curve, 6 = Oci^p), at the polar cap, where the argument varies from 
(p = — 7r/2 to (p = +71 /2. In the expected regime R^^/L <^ 1, the contribution Ax of the 
X-points to the solid angle of the polar cap is thus given by 



Ax-- 

Equation (IA2p implies 

cos[9c{(p)] 



+TT/2 j-ed-fi) 

dip / sin 9 d9 

tt/2 Jo 



-7r/2 



T 



n/2 



nl/2 



{1 - COS [9c{<p)]} dip 



(A3) 



cos (p 



R* 
2L 



cos If. 



(A4) 



Using this expression in equation ( IA3I) yields the approximation 



Ax^^. 



R* 

T 



(A5) 



From Figure HI we find that the other semi-curve defined by the magnetic field loops that 
connect one star to the other via z = +00 contributes (approximately) an equal amount. 
The solid angle of the polar cap thus becomes 



Ac 



2Rj, 



(A6) 



in agreement with the formula found empirically. 
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B. Generalized Electron Energy Distribution 

In this Appendix, we generalize the calculation of §2.4 to include a general power-law 
form for the energy distribution of fast electrons, i.e., 

n(7) d'j = nQ'y~^ d'j, (Bl) 

where the index p lies in the range 1 < p < 2.4. The lower limit leads to a flat spectral 
energy distribution uS^, ~ constant, whereas the upper limit is appropriate for collisionless 
shocks. 

The number density ut of fast particles, integrated from 7 = 1 to 7 = oo, is finite, and 
is given by ut = no/{p — 1) when p > 1. However, the integrated electron pressure is given 
by 

1 f°° 

Per = - / v-fmevn{-f) d-f, (B2) 



3 A 

which diverges unless p > 2. For p < 2, one must introduce an upper cuttoff in 7 to 
prevent the divergence of the pressure. For the case p > 2, with v ^ c, the pressure 
Pcr~nomecV[3(p-2)]. 



The volume emissivity of synchrotron radiation takes the generalized form 

-(p-l)/2 

where the dimensionless number Ci is defined by 



Cir.noeB^I—] , (B3) 



where T{Q is the Gamma function. For optically thin synchrotron emission, the flux density 
at a distance d from the source is given by 



Si, 



^ {CirenoeB±) L 



V 

v^) \d 



(B5) 



With 5„ oc zy-(p-i)/2 in equation (IBSp . the implied power Vy, integrated from zero frequency 
to z/, is given by 

Vy = ATcd^ r du = Aiid^ (B6) 
Jo V3 — P/ 

For the system parameters of V773 Tau A, we find the power Vu = (1.53 x 10^'^) /{3 —p) erg 
s~^ for synchrotron emission integrated from zero frequency to z/ = 90 GHz. Over a burst 
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time tB — 10^ s, the released energy £■ fa (1.53 x 10^°)/ (3 — p) at radio frequencies up to 90 
GHz. 

Using = 2 G, = (127)^ and = 400 mJy at 90 GHz, we find that the total 
number ut of fast electrons is given by 

nr^^^ -^^(127)(--%m-^ (B7) 

Table 1 presents the number density Ht and the parameter Ci for a range of indices p. Note 
that for 2 < p < 2.4, the density of fast electrons is consistent with the requirements arising 
from the reconnection-time argument (see text). 

The mean synchrotron optical depth through the source region at frequency u is given 

by 

1„ fcrA /^X-(P+4)/2 



where 



^ 3<-'./.r (!?L±1) r (^) . (B9, 

Table 1 also presents the optical depth t„ (at 1^ — 90 GHz) and the parameter C2 for varying 
indices p. For all plausible energy distributions of the fast electrons, those with p > 1, the 
synchrotron radiation is optically thin at the observed radio frequency z/ = 90 GHz. Note 
that radio observations at lower frequencies will provide important additional information 
regarding the physical conditions in the system (see text). 



Table 1. Electron Density and 90 GHz Optical Depth 



P Ci C2 riT [cm ^] T,, 



1.1 


7.94 


6.68 


1.92 X 10^ 


2.41 X 10" 


-4 


1.5 


5.52 


7.17 


3.83 X 10^ 


7.44 X 10- 


-5 


2.0 


4.42 


8.38 


2.70 X 10^ 


5.44 X 10" 


-5 


2.4 


4.12 


9.90 


1.44 X 10^ 


4.92 X 10- 


-5 
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